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ABSTRACT 

Pulmonary arterial hypertension (PAH) is a progressive and fatal disease of the lung vasculature for which the molecular etiologies 
are unclear. Specific metabolic alterations have been identified in animal models and in PAH patients, though existing data focus 
mainly on abnormalities of glucose homeostasis. We hypothesized that analysis of the entire metabolome in PAH would reveal 
multiple other metabolic changes relevant to disease pathogenesis and possible treatment. Layered transcriptomic and metabolomic 
analyses of human pulmonary microvascular endothelial cells (hPMVEC) expressing two different disease-causing mutations in 
the bone morphogenetic protein receptor type 2 (BMPR2) confirmed previously described increases in aerobic glycolysis but also 
uncovered significant upregulation of the pentose phosphate pathway, increases in nucleotide salvage and polyamine biosynthesis 
pathways, decreases in carnitine and fatty acid oxidation pathways, and major impairment of the tricarboxylic acid (TCA) cycle 
and failure of anaplerosis. As a proof of principle, we focused on the TCA cycle, predicting that isocitrate dehydrogenase (IDH) 
activity would be altered in PAH, and then demonstrating increased IDH activity not only in cultured hPMVEC expressing mutant 
BMPR2 but also in the serum of PAH patients. These results suggest that widespread metabolic changes are an important part 
of PAH pathogenesis, and that simultaneous identification and targeting of the multiple involved pathways may be a more fruitful 
therapeutic approach than targeting of any one individual pathway. 
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Pulmonary arterial hypertension (PAH) is a fatal, 
progressive disease of the pulmonary vasculature 
characterized by increasing pulmonary vascular resistance 
that leads to right heart failure and death. 1121 The disease 
exists in several forms in humans, including a heritable 
form caused primarily by mutations in bone morphogenetic 
protein receptor type 2 (BMPR2) and an idiopathic 
form that is clinically and in many ways molecularly 
indistinguishable from the inherited disease. 13 51 
Despite extensive investigations in PAH patients and 
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in a variety of animal models of PAH, the molecular 
mechanisms of disease pathogenesis have remained 
relatively obscure. Multiple converging lines of evidence 
point to disruption of interdependent metabolic pathways 
as being central to the molecular pathogenesis of PAH. 
In expression arrays from Bmpr2 mutant mice, nearly 
50% of the significantly altered genes fall into metabolic 
gene ontology groups, without identification of specific 
metabolic pathways. 161 Several animal models of PAH 
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show a shift toward aerobic glycolysis, the so-called 
"Warburg effect" that has been identified as central to 
malignant transformation in a number of tumor types. 17 " 91 
Alterations in glucose uptake and utilization, alongside 
changes in mitochondrial oxidative phosphorylation, 
have been demonstrated in the pulmonary artery 
endothelium from patients with PAH. 110 ' 111 More recently, 
PAH patients not previously known to have diabetes or 
any other obvious metabolic diseases were found to have 
measurable increases in hemoglobin Ale compared to 
age- and BMI-matched controls, suggesting that whole- 
body glucose homeostasis is impaired in PAH. 112,131 
Pulmonary hypertension associated with chronic hypoxia 
has been directly linked to an imbalance between 
glycolysis, glucose oxidation, and fatty acid oxidation. 191 
Finally, therapies aimed at normalizing glucose oxidation 
directly (e.g., inhibitors of pyruvate dehydrogenase 
kinase such as dichloroacetate] or via modulation of 
the balance between fatty acid oxidation and glucose 
oxidation (e.g., partial fatty acid oxidation inhibitors such 
as trimetazidine or ranolazine] have shown great promise 
in treating PAH and have demonstrated the importance 
of metabolic disturbances in disease initiation and 
maintenance. 18 ' 14181 Indeed, dichloroacetate has entered 
Phase I trials in humans (ClinicalTrials.gov identifier 
NCT01083524). Though the weight of evidence suggests 
that metabolic reprogramming is a key feature of the 
molecular pathogenesis of PAH, existing data focus mainly 
on abnormalities of glucose homeostasis, and the full 
breadth and scope of the altered metabolic pathways in 
PAH are unknown. 

We hypothesized that a broad-based metabolomic 
analysis of BMPR2 mutations that are known to cause 
PAH would reveal multiple coexisting and interdependent 
metabolic abnormalities beyond changes in glucose 
homeostasis. We quantify several hundred small molecule 
metabolites in native human pulmonary microvascular 
endothelial cells (hPMVEC] and in hPMVEC expressing 
one of two different disease-causing BMPR2 mutations. 
Organization of the significantly changed metabolites 
into known biochemical pathways confirms that multiple 
interconnected metabolic pathways are deranged in 
PAH. Gene expression array analysis from these same 
cells shows that metabolic genes represent the largest 
single group of significantly changed genes and support 
the findings from the metabolomic analyses. Using these 
layered metabolomic and transcriptomic analyses, we 
then predict alteration of the activity of a specific enzyme 
in the tricarboxylic acid (TCA] cycle - namely, isocitrate 
dehydrogenase (IDH] — as a proof of principle and 
demonstrate increased IDH activity in mutant hPMVEC 
and in the serum of patients with PAH. 



MATERIALS AND METHODS 

Human pulmonary microvascular endothelial 
cell culture 

Human PMVEC were grown in culture as previously 
described. 119 211 Cells were maintained in Endothelial Cell 
Growth Medium MV from PromoCell (Heidelberg, Germany] 
in standard cell culture incubators (37°C, humidified, 5% 
C0 2 ) and were used at or before the 10th passage. 

Generation of stably transfected hPMVEC 

Cells were transfected with either empty vector (native) 
or vector containing BMPR2 with R332X (KD] or 2579- 
2580delT (CD] mutations and stably selected using G418S 
as previously described. 1221 Endothelial character of the cells 
for this study was confirmed by immunohistochemistry for 
the von Willebrand factor and by analysis of expression 
arrays for a panel of endothelial markers (Figure S1A, B). 

Transcriptomic analysis 

Native and mutant hPMVEC were grown to 80% confluence, 
transitioned from G418S selection for at least 12 hours, and 
mRNA was isolated as described. 1231 Two AffymetrixHGU133 
Plus 2 arrays were run for each condition, with RNA for each 
array representing a pool of three independently grown 
plates, for a total of six arrays representing 18 biologically 
distinct events (three conditions x three plates x two arrays 
each). Results were analyzed using dChip and R statistical 
software. Significantly changed genes were determined 
using a requirement of a minimum of a 2x change, a 
minimum difference in expression of at least 200 arbitrary 
Affymetrix units, and a P<0.01 by a t-test for differences. 

Metabolomic analysis 

Full details of the methodology for the mass spectrometry- 
based metabolomic analyses are given in Supplemental 
Methods and as described previously 124 251 Briefly, samples 
(N=7 for each condition] were subjected to methanol 
extraction, split into aliquots for analysis by ultrahigh 
performance liquid chromatography/mass spectrometry 
(UHPLC/MS] in either the positive or negative ion mode 
or by gas chromatography/mass spectrometry (GC/MS]. 
Internal standards and controls for signal blank, technical 
replicates, and instrument performance were spiked into 
the samples and tracked throughout the analysis. Metabolite 
concentrations were determined by automated ion 
detection, manual visual curation, and were analyzed in-line 
using software developed by Metabolon. 1261 Significance was 
set at P<0.05 by Welch's two-sample t-test with correction 
for multiple comparisons using q-values. 1271 

NADP + -dependent IDH activity assays 

IDH activity assays were performed using the BioVision 
Isocitrate Dehydrogenase Activity Assay Kit (Mountain 
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View, Calif.) according to the manufacturer's instructions. 
For hPMVEC, cells were grown in 6-well plates to 
50 — 60% confluence to yield approximately 500,000 
cells per well, harvested, and lysed directly in the assay 
buffer. Aliquots were used for the IDH activity assay and 
for protein concentration determination by Pierce BCA 
assay. For serum, samples were used as undiluted 50 ul 
aliquots and assayed for IDH activity according to the 
instructions. 

Human subjects 

All patients and normal volunteers provided written 
informed consent to participate in research protocols 
approved by the institutional review boards of all 
participating institutions (IRB protocol number 9401). 
Blood was drawn by standard venipuncture, centrifuged 
to collect serum, and serum was stored at -80°C until 
analysis. 

Statistical analyses 

Analyses were performed using R statistical software and 
using GraphPad Prism. Welch's t-test or two-way ANOVA 
were used for tests of statistical significance. Box-and- 
whisker plots represent 25th — 75th percentiles with the 
box, the median with the center line, and Tukey whiskers 
representing 1.5 times the interquartile range. Scatter plots 
show individual data points with mean±SEM depicted. For 
most analyses, significance was set at P<0.05, with P<0.01 
being used as the significance threshold for RNA expression 
microarray analysis. 

RESULTS 

BMPR2 mutations resulted in widespread 
changes in endothelial cell gene expression that 
organized into specific pathways 

We sought to compare expression arrays from native 
hPMVEC to those from hPMVEC expressing one of 
two mutant BMPR2 constructs, and to organize the 
significantly different genes into functional pathways. 
Human pulmonary microvascular endothelial cells were 
stably transfected with either R332X mutation in the 
kinase domain (KD) or a 2579-2580delT mutation in the 
cytoplasmic tail domain (CD). The CD mutation has been 
previously shown to dysregulate BMPR2 interaction with 
and signaling through LIMK-1, c-Src, and Tctex-1; [28 311 the 
KD mutation also includes dysregulated signaling through 
the canonical Smad pathway. 1321 

Using a requirement of a minimum of twofold change, a 
minimum difference in expression of at least 200 arbitrary 
Affymetrix units, and aP<0.01 by t-test for difference, we 
found 687 probe sets representing 507 unique Entrez IDs, 
with common changes between both BMPR2 mutants 



and native hPMVEC, with a false discovery rate (FDR) 
of zero (determined by scrambling group identifiers). 
These data have been deposited in GEO, accession number 
pending , and a full list of the 507 genes is provided in 
Supplemental Dataset SI. Distribution of gene ontology 
groups was nearly identical to our previously published 
expression arrays interrogating PMVEC isolated from 
our Bmpr2 R899x and Bmpr2 delx4+ mouse models. 16 33341 
These included genes involved in apoptosis, proliferation, 
stimulus response, cytoskeletal organization, and 
development (Fig. 1). Roughly 40% of the genes changed 
(216/507) were related to small molecule metabolism. 
Heterologous expression of BMPR2 mutations resulted 
in broad changes in TCA cycle, glycolysis, hypoxia- 
inducible factor (HIF) responsive metabolic elements, and 
carnitine, fatty acid, and glutamate metabolism compared 
to expression of native BMPR2. Relatively unaffected 
pathways include glycan synthesis and metabolism, 
vitamin/cofactor metabolism (with the exception of 
folate and single-carbon metabolism), and xenobiotic 
metabolism. Thus, the affected pathways showed a degree 
of specificity as opposed to nonspecific whole metabolome 
dysfunction. 

Metabolomic analysis of BMPR2 mutant 
endothelial cells showed significant alteration 
of multiple interdependent metabolic pathways 

To determine the whole metabolome consequences 
of disease-causing BMPR2 mutation in endothelial 
cells, we undertook a simultaneous multiplexed mass 
spectrometric quantification of several hundred small 
molecule metabolites in CD and KD mutant hPMVEC 
and compared these mutations to the native hPMVEC. 
In this analysis, 267 small molecule metabolites were 
confidently identified in seven biological replicates for 
each condition described above (native, CD, and KD, 
Figure S2). Significantly changed biochemicals from the 
native condition were identified as those biochemicals 
with a P-value <0.05 based upon Welch's two-sample 
t-test, which had a maximum FDR of 3.2% based upon 
q-values [271 for that set of biochemicals withP-values <0.05. 
The full dataset is provided in Supplemental Dataset S2. 
Compared to the native hPMVEC, the CD mutants showed 
significant changes in 65% of the metabolites quantified 
(172/267, with 87 increased and 85 decreased) and the 
KD mutants showed significant changes in 37% of the 
metabolites (99/267, with 61 increased and 38 decreased). 
This represented confident identification of approximately 
11% of the database of named compounds available in 
this analysis, with the CD mutants showing significant 
changes in the levels of approximately 7% of the total 
compounds in the database. For the KD mutants, a further 
14% (38 metabolites) approached statistical significance 
(0.05<P<0.10 by Welch's two-sample t-test), though these 
were not included in the analyses discussed below. 
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Figure 1 : Transcriptomic analysis of human pulmonary 
microvascular endothelial cells expressing one of two 
dominant negative BMPR2 mutant constructs. The pie 
chart depicts the pathways represented by the 507 genes 
showing significantly altered expression in the CD and 
KD mutant hPMVEC compared to the native hPMVEC. 
The largest single major ontology group was metabolic 
genes, and the heatmap at the right of the figure shows 
changes in individual genes within this larger group 
broken down by subsets, with red indicating increased 
expression and blue denoting decreased expression. 



Increased pentose phosphate pathway 
metabolites and polyamine biosynthesis indicate 
increased proliferation in BMPR2 mutant 
endothelial cells 

Intermediates in the pentose phosphate pathway for BMPR2 
mutant hPMVEC along with changes in the expression of the 
corresponding genes compared to native hPMVEC are shown 
(Fig. 2]. The pentose phosphate pathway interfaces with 
multiple other metabolic pathways, including the glycolytic 
pathway and NADPH synthetic pathways, in addition to 
providing 5-carbon sugars for nucleotide synthesis. The 
enzyme primarily responsible for synthesis of NADPH 
in the pentose phosphate pathway, glucose-6-phosphate 
dehydrogenase, exhibited significantly reduced expression 
in both the CD and KD mutants. The upregulation of the 
pentose phosphate pathway in the BMPR2 mutant hPMVEC 
thus was apparently driven by upregulation from ribose-5- 
phosphate isomerase downstream that overcame decreases 
in glucose-6-phosphate dehydrogenase expression and 
that directly related to increased nucleotide synthesis and 
salvage. The significantly increased levels of purine and 
pyrimidine nucleosides downstream from both mutations 
further supported this conclusion. Finally, intermediates in 
the terminal portion of the polyamine synthesis pathway in 
the CD and KD mutants were increased preferentially over 
urea cycle intermediates (Figure S3], further supporting the 



increased cell proliferation previously observed for these 
specific mutations. 161 

Multiple energy metabolism pathways were 
disrupted in BMPR2 mutant endothelium 

Several groups have previously shown an increase in 
aerobic glycolysis (the "Warburg effect"] in PAH, [7 9,35] and 
increased glucose uptake has been demonstrated in the 
lungs of PAH patients by positron emission tomography. 1101 
In this whole metabolome analysis, glycolysis showed 
significant upregulation as a consequence of BMPR2 
mutation, though this was not to the same extent for 
both mutations (Fig. 3). Glycolytic metabolites were 
significantly increased in the CD mutants for all of the 
pathway intermediates down to pyruvate, whereas the 
KD mutants showed statistically significant increases or 
trends in increases only for the metabolites earlier in the 
glycolytic pathway. Consistent with global activation of the 
glycolytic pathway, expression of the majority of the genes 
for glycolytic enzymes was increased, including genes 
coding for glucose uptake transporters. The genes for two 
of the three regulated enzymes in the pathway, hexokinase, 
and phosphofructokinase, showed decreased expression in 
arrays from the BMPR2 mutant hPMVEC, though the third 
regulated enzyme, pyruvate kinase, showed increased 
expression for at least one isoform. 
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Figure 2: Intermediates in the pentose phosphate pathway specific to nucleotide synthesis are increased in BMPR2 mutant hPMVEC. Major intermediates in 
the pentose phosphate pathway are shown. In all graphs, native hPMVEC are in white boxes, CD hPMVEC in vertical hatched boxes, and KD hPMVEC in 
diagonal hatched boxes. Quantities are in arbitrary units specific to the internal standards for each quantified metabolite and normalized to protein concentration. 
N = 7 for each box, with whiskers indicating Tukey whiskers and extreme data points indicated by filled circles. *P<0.05 compared to native. Genes coding 
for the enzymes that catalyze particular steps in the pathway are indicated by their Entrez Gene names, with red indicating significantly increased expression 
in the transcriptomic analysis and blue indicating significantly decreased expression. The four graphs at the bottom of the figure show quantitation of purine 
and pyrimidine nucleosides. G6PD, glucose-6-phosphate dehydrogenase; H6PD, hexose-6-phosphate dehydrogenase; RPIA, ribose-5-phosphate isomerase A; 
RPE, ribulose-5-phosphate-3-epimerase; TKT, transketolase. 



The interplay between fatty acid oxidation and glucose 
utilization has been shown to play an important role in 
pulmonary hypertension related to chronic hypoxia and 
in right ventricular hypertrophy and failure induced by 
pressure overload in a pulmonary artery banding model, 19181 
but this has not been explored in pulmonary arterial 
hypertension specifically. We thus sought evidence for 
alterations in the major pathways for fatty acid oxidative 
metabolism in the context of disease-causing BMPR2 
mutations. We found that carnitine and its downstream acyl 
metabolites were significantly reduced in the CD and KD 
mutant hPMVEC compared to the native endothelial cells 
(Fig. 4). Decreased levels of carnitine itself as well as glycine 
(a by-product of carnitine synthesis] suggested decreased 
synthesis of carnitine itself. Levels of palmitoylcarnitine, 



isobutyrylcarnitine, and propionylcarnitine were also 
significantly decreased. Decreased expression of many of the 
key genes involved in carnitine/acylcarnitine metabolism 
and trafficking was also observed, including the two 
major carnitine palmitoyltransferase genes and one of the 
major carnitine/acylcarnitine translocases. We also found 
significantly decreased expression of a number of the acyl- 
CoA dehydrogenase genes involved in fatty acid oxidation. 

Activity of the tricarboxylic acid [TCA] cycle has been 
shown to be reduced in a variety of different types of cancer, 
and this has been proposed to be a central advantage 
exploited by cancer cells, allowing for diversion of TCA 
cycle intermediates toward macromolecule synthesis 
while relying on other energy-generating pathways such as 
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Figure 3: Glycolysis is significantly upregulated in BMPR2 mutant hPMVEC, particularly in cytoplasmic tail domain mutants. The classical glycolysis 
pathway intermediates are shown. In all graphs, native hPMVEC are in white boxes, CD hPMVEC in vertical hatched boxes, and KD hPMVEC in diagonal 
hatched boxes. Quantities are in arbitrary units specific to the internal standards for each quantified metabolite and normalized to protein concentration. N = 7 
for each box, with whiskers indicating Tukey whiskers and extreme data points indicated by filled circles. */ > <0.05 compared to native. Genes coding for the 
enzymes that catalyze particular steps in the pathway are indicated by their Entrez Gene names, with red indicating significantly increased expression in the 
transcriptomic analysis and blue indicating significantly decreased expression. SLC2A1 and SLC2A3, solute carrier family 2 (facilitated glucose transporter), 
members 1 and 3; HK2, hexokinase 2; PFKL and PFKM, phosphofructokinase, liver and muscle isoforms; ALDOA, aldolase A; TPL1, triosephosphate isomerase 
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glycolysis. 136 381 Although a number of the metabolic features 
of cancer cells have been observed in PAH, defects in the 
TCA cycle have not been extensively described. In hPMVEC 
expressing BMPR2 mutations, there were extensive 
metabolic defects in the TCA cycle indicating overall 
decreased activity of the cycle downstream from citrate 
(Fig. 5). In particular, significant decreases in succinate, 
fumarate, and malate were present in the CD mutants, 
whereas much more mild nonsignificant decreases in the 
mean concentrations of succinate and malate were observed 
for the KD mutants. In both mutants, concentrations of 
citrate were significantly increased, and concentrations of 
pyruvate and lactate were equivalent to the native hPMVEC, 
suggesting that the defect in the TCA cycle occurred distal to 



citrate. Moreover, this suggested that alternative catabolic 
pathways (e.g., fatty acid oxidation, peptide/amino acid 
catabolism] feeding into the TCA cycle were insufficient 
to support concentrations of succinate, fumarate, and 
malate, particularly in the CD mutants. The balance of 
TCA cycle intermediates is normally maintained by the 
complementary processes of anaplerosis and cataplerosis. 
Broadly defined, anaplerosis refers to the addition of 
4- and 5-carbon intermediates into the TCA cycle (e.g., 
oxaloacetate, alpha-ketoglutarate, and succinyl-CoA] using 
pyruvate, aspartate, glutamate, or fatty acids as substrates, 
to support mitochondrial respiration and to replenish TCA 
cycle intermediates diverted to biosynthesis. Cataplerosis 
refers to the removal of these same TCA cycle intermediates 
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Figure 4: Carnitine metabolism and fatty acid oxidation are significantly depressed in BMPR2 mutant hPMVEC. Multiple carnitine metabolites and their flow 
into fatty acid oxidation are shown. Intermediates for which significant differences in one or both mutant conditions were detected are shown. In all graphs, 
native hPMVEC are in white boxes, CD hPMVEC in vertical hatched boxes, and KD hPMVEC in diagonal hatched boxes. Quantities are in arbitrary units 
specific to the internal standards for each quantified metabolite and normalized to protein concentration. N = 7 for each box, with whiskers indicating Tukey 
whiskers and extreme data points indicated by filled circles. *P<0.05 compared to native. Genes coding for the enzymes that catalyze particular steps in the 
pathway are indicated by their Entrez Gene names, with red indicating significantly increased expression in the transcriptomic analysis and blue indicating 
significantly decreased expression. CPT1A and CPT2, carnitine palmitoyltransferase isoforms 1A and 2; SLC25A20, camitine/acylcarnitine translocase; 
ACADS, ACADM, ACADSB, and ACADVL, acyl-CoA dehydrogenases - short chain, medium chain, short/branched chain, and very long chain; MLYCD, 
malonyl-CoA decarboxylase. 



to support biosynthetic processes such as gluconeogenesis, 
glyceroneogenesis, and fatty acid synthesis. [39] 

Two key specific anaplerotic pathways are the conversion 
of glutamine to glutamate and then to alpha-ketoglutarate, 
and the conversion of aspartate to oxaloacetate. [40] 
Levels of these specific amino acids were quantified 
and found to be significantly reduced in the CD mutant 
hPMVEC compared to native cells (Fig. 6). By contrast, 
the KD mutants showed increased levels of glutamine and 
glutamate, which likely contributed to the more modest 
reductions of TCA cycle intermediates in these cells, as 
there was more glutamine and glutamate available for 
the anaplerotic synthesis of alpha-ketoglutarate. The 



differences between the CD and KD mutants may be at 
least in part attributable to differences in the synthesis 
of N-acetylaspartylglutamate (NAAG). The CD mutants 
showed increased expression of NAAG synthetase 
and increased concentrations of NAAG, whereas the 
KD mutants showed decreased expression of NAAG 
synthetase and a subsequently lower level of NAAG 
compared to the CD mutants. It thus is likely that, in the CD 
mutants, much of the glutamate and aspartate that might 
otherwise have been used to feed into the TCA cycle was 
instead being used for the synthesis of NAAG. This was 
not the case in the KD mutants. Both mutations appeared 
to drive an overall increase in peptide and amino acid 
catabolism, as both mutants showed significant increases 
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Figure 5: TCA cycle intermediates are significantly decreased in BMPR2 mutant hPMVEC, particularly in CD mutants. The TCA cycle and its major 
intermediates are shown. In all graphs, native hPMVEC are in white boxes, CD hPMVEC in vertical hatched boxes, and KD hPMVEC in diagonal hatched 
boxes. Quantities are in arbitrary units specific to the internal standards for each quantified metabolite and normalized to protein concentration. N = 7 for each 
box, with whiskers indicating Tukey whiskers and extreme data points indicated by filled circles. *P< 0.05 compared to native. Genes coding for the enzymes 
that catalyze particular steps in the pathway are indicated by their Entrez Gene names, with red indicating significantly increased expression in the transcriptomic 
analysis and blue indicating significantly decreased expression. PDK1-4, pyruvate dehydrogenase kinase 1-4; PDHA1, pyruvate dehydrogenase (lipoamide) 
alpha 1 ; PDHB, pyruvate dehydrogenase El component, beta subunit; DLAT, dihydrolipoamide S-acetyltransferase; AC02, aconitase 2; IDH3A/B/G, isocitrate 
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in the concentrations of dipeptides (Figure S4], though 
this was more pronounced in the CD mutants. In addition, 
anaplerosis via branched chain amino acid metabolism 
appeared to be more significantly impaired in the CD 
mutants compared to the KD mutants, as evidenced 
by more significant decreases in isobutyrylcarnitine 
and propionylcarnitine, metabolites that participate in 
branched chain amino acid metabolism as well as fatty 
acid metabolism. The CD mutants exhibited a major 
failure of anaplerosis on multiple levels that was much 
more mild in the KD mutants, and thus must rely on 
multiple "salvage" pathways of peptide/amino acid 
catabolism. 



Predicted differences in IDHT/2 activity were 
present in mutant hPMVEC and in human PAH 
patients 

We chose to examine the TCA cycle in more detail, as 
this is a major point of integration of multiple pathways 
involved in energy production as well as biosynthesis and 
so might better reflect the summation of alterations in 
these many different pathways. On closer inspection of the 
transcriptomic and metabolomic data relevant to the TCA 
cycle, there was a clear functional change in the BMPR2 
mutants compared to the wild-type endothelial cells that 
occurred somewhere between citrate and succinate. The 
two most likely enzymatic candidates were aconitase and 
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Figure 6: Glutamine/glutamate and aspartate metabolism, two major anaplerotic pathways, are significantly reduced in CD mutant hPMVEC. Intermediates 
for which significant differences in one or both mutant conditions were detected are shown. In all graphs, native hPMVEC are in white boxes, CD hPMVEC 
in vertical hatched boxes, and KD hPMVEC in diagonal hatched boxes. Quantities are in arbitrary units specific to the internal standards for each quantified 
metabolite and normalized to protein concentration. N = 7 for each box, with whiskers indicating Tukey whiskers and extreme data points indicated by filled 
circles. *P<0.05 compared to native. Genes coding for the enzymes that catalyze particular steps in the pathway are indicated by their Entrez Gene names, with 
red indicating significantly increased expression in the transcriptomic analysis and blue indicating significantly decreased expression. GLUD1 and GLUD2, 
glutamate dehydrogenase 1 and 2; GLS, glutaminase; GLUL, glutamate-ammonia ligase; RIMKLB, N-acetylaspartylglutamate synthetase B. 



isocitrate dehydrogenase (IDH). Of these two enzymes, 
aconitase has been shown to be inactivated by oxidative 
stress, 1411 and oxidative stress is known to be increased 
in the context of BMPR2 mutations. [2Z] To remove this as 
a confounding factor, we chose to examine IDH activity. 
While NAD + -dependent IDH activity (corresponding to 
the IDH3 isoform] did not differ between wild-type and 
BMPR2 mutant endothelial cells, NADP + -dependent IDH 
activity (corresponding to IDH isoforms 1 and 2, hereafter 
IDH 1/2) was significantly increased in BMPR2 mutant 
endothelial cells (Fig. 7A). We then sought to determine 
if these findings were applicable to patients with PAH. 
We quantified IDH1/2 activity in serum from controls, 
from patients with heritable PAH known to have BMPR2 
mutations, and from two pooled cohorts of patients with 
IPAH (one from the United States and one from Japan). 



Serum IDH 1/2 activity was significantly increased in both 
HPAH and IPAH patients compared to controls (Fig. 7B). 
The variability in IDH1/2 activity observed in the PAH 
patients likely does reflect variability in disease activity, at 
least in part. We separated the PAH patients for whom data 
were available into two groups based upon the presence 
or absence of any prostanoid therapy (intravenous, 
subcutaneous, or inhaled prostacyclin receptor agonist) 
and analyzed serum IDH1/2 activity (Figure S5). Though 
only approaching statistical significance (P=0.1 by Welch's 
t-test), the mean serum IDH 1/2 activity for the prostanoid 
treated group of PAH patients clearly trended toward 
being higher than for the nontreated group. As prostanoid 
therapy is only initiated in more severe disease, the higher 
serum IDH activity may well be a reflection of disease 
severity. 
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Figure 7: IDH1/2 activity is increased in mutant hPMVEC and in serum from patients with pulmonary arterial hypertension. (A) Compared to native hPMVEC 
(1.07±0.04 mU/5><10 5 cells), the CD and KD mutant-expressing hPMVEC both show significantly increased NADP~-dependent IDH activity in cell lysates 
(1.56±0.04 and 1.38 +/0.07 mU/5xl0 5 cells, respectively). N=6 for each group, mean±SEM indicated, /><0.0001 by one-way ANOVA, **P<0.0001 vs. native, 
*P<0.003 vs. native. (B) NADP + -dependent IDH activity in the serum from patients with heritable PAH (N=8, 1.83±0.1 mU/mL) and idiopathic PAH (N=17, 
1.85±0.09 mU/mL) was significantly increased compared to serum from normal control individuals (N=13, 1.61±0.04mU/mL), mean±SEM indicated, **P<0.03 
vs. control by Welch's t-test. 



DISCUSSION 

We present here a whole metabolome analysis of the 
effects of BMPR2 mutations known to cause pulmonary 
arterial hypertension. We have constructed an integrated 
picture of the complex and interdependent metabolic 
changes that occur downstream from BMPR2 mutations 
in human pulmonary microvascular endothelial cells, and 
this integrated view has demonstrated more widespread 
metabolic defects in PAH than have been previously 
known. The shift toward aerobic glycolysis that is typified 
by the Warburg effect in PAH has been known for some 
time. 171 More recently the interplay between fatty acid 
oxidation and control of glycolysis versus glucose oxidation 
(typified by the Randle cycle) has been demonstrated to 
be important in hypoxic pulmonary hypertension. [942] 
In addition to confirming these previously identified 
metabolic defects in PAH, we have identified significant 
alterations in the TCA cycle and at least some of the 
pathways that interact with it. 

Our analysis reveals a profound failure of anaplerosis 
present downstream from BMPR2 mutations that has been 
largely unexplored as a mechanism of disease pathogenesis 
in PAH. TCA cycle intermediates are depleted through what 
appears to be a combination of decreased activity of the 
cycle itself plus abnormal shunting of intermediates from 
other pathways (e.g., glutamate, glutamine, aspartate, and 
branched chain amino acids] that could otherwise be used 
to replenish the intermediates of the TCA cycle. This implies 
that the metabolic defects in PAH cannot be simply summed 
up under the umbrella of Warburg metabolism. [3543] It 
is possible that diversion of TCA cycle intermediates for 
biosynthesis and reliance on aerobic glycolysis is actually 



a feature of decreased BMP signaling, as would be seen 
in settings of tissue remodeling and repair, for example. 
However, these would be settings of temporary loss or 
reduction of BMP signaling, and what might be adaptive 
cataplerosis in the short-term becomes pathogenic in the 
setting of permanently decreased BMP signaling due to 
mutation. More importantly, the therapeutic implication is 
that it is very likely that most or all of the errant pathways 
will need to be targeted to expect a significant therapeutic 
impact, and current investigational therapies address only 
certain dysfunctional pathways that largely do not address 
the failure of anaplerosis. 

To demonstrate further the utility of our combined 
metabolomic and transcriptomic analysis, we chose to 
quantify NADP + -dependent IDH activity. We predicted 
that this enzymatic activity would be altered in PAH and 
then demonstrated the accuracy of that prediction in both 
cultured cells and in patients with PAH. Alterations in 
IDH activity have not been previously described in PAH; 
however, multiple conceptual links between PAH and 
cancer have been proposed, and there is a growing body of 
literature linking altered IDH activity in a causative way to 
at least certain types of cancer, [2444 " 47] so the identification 
of altered IDH1/2 activity in PAH may in fact be directly 
related to disease pathogenesis. 

A second possibility is that the anaplerotic failure and 
increased IDH 1/2 activities in PAH are driving forces 
for HIF activation that underlies PAH, [1548] which then 
perpetuates increased aerobic glycolysis, apoptosis 
resistance, and decreased mitochondrial number. [11] HIF 
can be activated by decreased alpha-ketoglutarate and 
increased citrate concentrations, both of which activate 
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HIF by decreasing the efficiency of prolyl hydroxylase. [49] 
Under hypoxic conditions, activation of HIF has very 
recently been shown to increase IDH2-mediated 
conversion of glutamine-derived alpha-ketoglutarate to 
citrate. [50] The resultant hypothesized increase in citrate 
and decrease in alpha-ketoglutarate brought about by 
increased IDH2 activity would thus be predicted to further 
drive HIF activation, setting up something of a vicious 
cycle. Our data show significant changes in the expression 
of HIF responsive genes in BMPR2 mutant hPMVEC. 
Alternatively, HIF activation can be driven by increased 
oxidative stress and decreased antioxidant defenses, as 
is seen with epigenetic inactivation of SOD2, [511 and this 
may be the upstream event that drives IDH activation. HIF 
activation can drive metabolic reprogramming that leads 
to increased IDH activity, [50] and IDH may be upregulated 
in response to oxidative stress to serve as a source of 
NADPH that is used to maintain reduced glutathione pools 
intracellularly. [52] 

More importantly, IDH activity appears to track with 
disease activity, as patients with more severe disease (i.e., 
those treated with a prostacyclin agonist] trend toward 
higher serum IDH activity. The fact that those patients with 
disease severe enough to warrant prostanoid therapy still 
exhibit increased serum IDH activity after the initiation 
of therapy suggests that our most efficacious class of 
drugs for PAH still does not correct all of the underlying 
metabolic defects in PAH and highlights the need for 
therapies that arise from a deeper understanding of the 
molecular pathogenesis of PAH. Further investigations 
are needed to determine if the increased IDH activity 
described in PAH in this study is truly pathogenic, adaptive, 
or epiphenomenon. 

This type of investigation is not without limitations. The 
pulmonary endothelial cell clearly plays an important role 
in the pathogenesis of PAH, but many other cell types, 
including smooth muscle cells, 1531 lung mesenchymal stem 
cells, 1541 and resident immune cells 1551 all likely contribute 
to disease development. It is not possible in most cases in 
this study to determine what changes are truly causative 
of disease and which are consequences of the disease, 
although examination of cells in culture minimizes 
this problem to some degree. [23] While identification 
of putative molecular targets is possible and is greatly 
enhanced by this type of layered analysis, each target 
must be further validated independently and placed into 
a context. Still, this study has permitted the identification 
of previously unrecognized pathways that likely directly 
contribute to the development of PAH; the identification 
of promising new biomarkers, of which IDH1/2 activity 
is but one, to guide diagnosis and therapeutic evaluation; 
and the recognition of the importance of defining and 
simultaneously targeting the multiple affected metabolic 



pathways in future therapeutic development. We 
hypothesize that this approach would be similarly fruitful 
in many other complex diseases. 

APPENDIX 

Supplemental Methods 

Metabolomic Analysis: Native and mutant hPMVEC were 
grown to 80% confluence (yielding approximately 1-3 
xlO 7 cells per sample], transitioned from G418 sulfate 
selection to complete media without antibiotic for at least 
12 hours, and washed with phosphate-buffered saline 
and trypsinized to harvest cells. Cells were pelleted, the 
supernatant removed, and the dry pellet snap frozen in 
liquid nitrogen and stored at -80C until analysis. The 
non-targeted metabolic profiling platform employed for 
this analysis combined three independent platforms: 
ultrahigh performance liquid chromatography/tandem 
mass spectrometry (UHPLC/MS/MS 2 ] optimized for basic 
species, UHPLC/MS/MS 2 optimized for acidic species, 
and gas chromatography/mass spectrometry (GC/ 
MS]. Samples were processed essentially as described 
previously. [12] For each sample, 100L was used for 
analyses. Using an automated liquid handler (Hamilton 
LabStar, Salt Lake City, UT], protein was precipitated 
from the homogenate with methanol that contained four 
standards to report on extraction efficiency. The resulting 
supernatant was split into equal aliquots for analysis 
on the three platforms. Aliquots, dried under nitrogen 
and vacuum-desiccated, were subsequently either 
reconstituted in 50L 0.1% formic acid in water (acidic 
conditions] or in 50L 6.5mM ammonium bicarbonate 
in water, pH 8 (basic conditions] for the two UHPLC/ 
MS/MS 2 analyses or derivatized to a final volume of 50L 
for GC/MS analysis using equal parts bistrimethyl-silyl- 
trifluoroacetamide and solvent mixture acetonitrile:dich 
loromethane:cyclohexane (5:4:1] with 5% triethylamine 
at 60°C for one hour. In addition, three types of controls 
were analyzed in concert with the experimental samples: 
aliquots of a well-characterized human plasma pool 
served as technical replicates throughout the data set, 
extracted water samples served as process blanks, 
and a cocktail of standards spiked into every analyzed 
sample allowed instrument performance monitoring. 
Experimental samples and controls were randomized 
across platform run days. 

For UHLC/MS/MS 2 analysis, aliquots were separated 
using a Waters Acquity UPLC (Waters, Millford, MA] 
and analyzed using an LTQ mass spectrometer (Thermo 
Fisher Scientific, Inc., Waltham, MA] which consisted of 
an electrospray ionization (ESI] source and linear ion- 
trap (LIT] mass analyzer. The MS instrument scanned 
99-1000 m/z and alternated between MS and MS 2 scans 
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using dynamic exclusion with approximately 6 scans per 
second. Derivatized samples for GC/MS were separated 
on a 5% phenyldimethyl silicone column with helium 
as the carrier gas and a temperature ramp from 60°C to 
340°C and then analyzed on a Thermo-Finnigan Trace 
DSQ MS (Thermo Fisher Scientific, Inc.] operated at unit 
mass resolving power with electron impact ionization 
and a 50-750 atomic mass unit scan range. 

Metabolites were identified by automated comparison 
of the ion features in the experimental samples to a 
reference library of chemical standard entries that 
included retention time, molecular weight [m/z], 
preferred adducts, and in-source fragments as well 
as associated MS spectra, and were curated by visual 
inspection for quality control using software developed 
at Metabolon. 131 

For statistical analyses and data display purposes, any 
missing values were assumed to be below the limits 
of detection and these values were imputed with the 
compound minimum (minimum value imputation]. 
Statistical analysis of log-transformed data was performed 
using "R" [http://cran.r-project.org/]. Welch's t-tests 
were performed to compare data between experimental 
groups. A p-value of < 0.05 was considered statistically 
significant and multiple comparisons were accounted 
for by estimating the false discovery rate [FDR] using 
q-values. [4] 
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Figure SI : Native and mutant hPMVEC retain endothelial character in culture. (A) Immunohistochemistry for von Willebrand factor (red, with DAPI nuclear 
stain in blue) shows that hPMVEC retain functional expression of a canonical endothelial protein, and that expression of the mutant BMPR2 constructs does 
not change this substantially. Cells shown are representative. (B) A panel of ten endothelial markers quantified in expression arrays shows that expression of 
these markers is retained in the mutant hPMVEC, confirming endothelial character. Markers are as follows: FLT1, vascular endothelial growth factor receptor 
i,ESMl, endothelial cell-specific molecule \ \MMRN2, multimerin 2; SELP, selectin P; KDR, vascular endothelial growth factor receptor 2; PECAM1, platelet/ 
endothelial cell adhesion molecule; THBD, thrombomodulin; ICAM2, intercellular adhesion molecule 2; ENG, endoglin; VEGF, vascular endothelial growth factor. 
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Figure S2: Overall heatmap for the metabolomic changes in CD and KD 
mutant hPMVEC compared to native hPMVEC. Metabolites are shown 
grouped by superpathway and subpathway. Data are displayed as fold change in 
each metabolite for each mutant relative to the native. Statistically significant 
changes are denoted with red for increase above the level in native hPMVEC 
and green for decrease below native levels. 
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Figure S3: The polyamine biosynthetic pathway is significant upregulated in BMPR2 mutant hPMVEC. The major intermediates in the urea cycle and the 
polyamine biosynthetic pathway are listed. Intermediates that were confidently identified at detectable quantities are shown. In all graphs, native hPMVEC are 
in white boxes, CD hPMVEC in vertical hatched boxes, and KD hPMVEC in diagonal hatched boxes. Quantities are in arbitrary units specific to the internal 
standards for each quantified metabolite and normalized to protein concentration. N=7 for each box, with whiskers indicating Tukey whiskers and extreme 
data points indicated by filled circles. *i 3 <0.05 compared to native. Genes coding for the enzymes that catalyze particular steps in the pathway are indicated 
by their Entrez Gene names, with red indicating significantly increased expression in the transcriptomic analysis and blue indicating significantly decreased. 
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Figure S4: Dipeptides are significantly increased in mutant hPMVEC. Data 
are displayed as fold change in each metabolite for each mutant relative to 
the native. Statistically significant changes are denoted with red for increase 
above the level in native hPMVEC. 
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Figure S5: Serum IDH1/2 activity trends toward being higher in PAH patients 
on prostanoid therapy compared to PAH patients not treated with a prostacyclin 
receptor agonist. IDH activity measurements from heritable and idiopathic 
PAH patients were pooled and grouped based on prostanoid therapy status 
at the time of sample collection for those patients for whom such data were 
available. N=13 in the prostanoid group, N=6 in the no prostanoid group, 
mean±SEM for each group indicated, P=0A by Welch's two-sample t-test. 
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